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Abstract. Novel methods for diagnostics of molecular hydrogen plasma processes, 
such as ionization, production of high vibrational levels, dissociation of molecules 
via excitation to singlet and triplet states and production of metastable states, are 
presented for molecular hydrogen plasmas in corona equilibrium. The methods are 
based on comparison of rate coefficients of plasma processes and optical emission 
spectroscopy of lowest singlet and triplet transitions, i.e. Lyman-band 

and molecular continuum ^ 6^E+), of the hydrogen molecule in 

VUV wavelength range. Comparison of rate coefficients of spin-allowed and/or 
spin-forbidden excitations reduces the uncertainty caused by the non-equilibrium 
distributions of electron energy and molecular vibrational level, which are typically 
known poorly in plasma sources. The described methods are applied to estimate the 
rates of various plasma processes in a filament arc discharge. 
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1. Introduction 

Molecular hydrogen plasmas play an important role in negative (H“) [1] and molecular 
(H^) hydrogen [2H1] ion sources. Furthermore, hydrogen molecules are an important 
species in magnetically conhned plasmas of fusion devices |5] as well as in astrophysics 

0 . 

Optical emission spectroscopy (OES) is a powerful and widely used plasma 
diagnostic method. Wavelength range from near ultraviolet to near infrared is often 
used for the purpose (see e.g. and references therein). Most of the diagnostic methods 
based on OES are developed to obtain information about plasma parameters, such as 
the plasma density and temperature, which requires employing sophisticated collisional 
radiative models. Very few methods based on molecular emission have been developed 
for diagnostics of hydrogen plasma processes. For example, Sawada and Fujimoto |H] 
have developed methods for determining ionization and dissociation rate coefficients 
of hydrogen molecules based on optical emission spectroscopy of hydrogen atoms in 
visible light range. Graham has estimated the production rate of ground state hydrogen 
molecules at high vibrational levels based on (absolute) Lyman-alpha emission and 
comparison of the rate coefficients of molecular (singlet) and atomic excitations [9]. 
Lavrov et ah have demonstrated that molecular continuum emission could be used as a 
probe for dissociation rate in molecular hydrogen plasma [lO]. However, only excitations 
to and states were taken into account and the total calculated emission based 

y ^ 

on an extrapolation of spectral measurement data in the range of 225-450 nm. 

Measurement of the light emission in the VUV wavelength range (100-250 nm) is a 
relevant optical diagnostic of certain plasma processes in molecular hydrogen plasmas. 
This is because the transitions from the lowest excited electronic states to the ground 
state emit photons in this range, radiative lifetimes of the lowest excited states are short 
and electron impact excitation cross sections from the ground state to those states are 
approximately an order of magnitude greater than excitation cross sections to upper 
electronic states. Therefore, the physical origin of the population density of the lowest 
excited states can be understood straightforwardly and photon emission rate is linearly 
proportional to the electron impact excitation rate from the ground state. 

Application of VUV-emission for diagnostics of molecular hydrogen plasma 
processes such as ionization, dissociation and excitation, is demonstrated in this paper. 
Instead of using information about the plasma parameters i.e. density and temperature, 
the applied method is based on direct measurement of the VUV-emission rate, which is 
proportional to the electron impact excitation rate. Comparison of the rate coefficients 
results to a robust plasma diagnostic method requiring minimum knowledge on the 
plasma parameters. The method can be applied to obtain absolute numbers or relative 
changes of volumetric production rates through absolute or relative measurements of 
the VUV-irradiance. The error caused by the electron energy distribution (EED) and 
molecular vibrational temperature is discussed. 

The presented methods are applied for a set of previously presented VUV-irradiance 


3 


measurements of a filament driven negative hydrogen ion source, LIISA, at the JYFL 
accelerator laboratory. Estimated excitation and ionization rates based on saturation 
value (at >100 eV electron energy) of rate coefficients of vibrationally cold plasma, 
have been presented and discussed briefly in the earlier study. In this study the analysis 
is carefully extended to account for different EEDF’s and vibrational temperatures. 
Furthermore, the dissociation rate of hydrogen molecules and the production rate of 
metastable hydrogen molecules (c^II^j state) are determined for the hrst time. 

2. Processes affecting the population densities of excited electronic states 
of hydrogen molecules 

Hydrogen molecule has two multiplet systems of electronic states, singlet and triplet, 
which differ by the orientation of the electron spins (iS'=0 for singlet and 5'=! for triplet) 
governed by the Pauli exclusion principle. Radiative transitions between the electronic 
states are classihed as optically allowed and optically forbidden. Optically allowed 
electron transitions are subject to the following selection rules: (a) change of the parity 
of the total wave function {g u), (b) conservation of the total spin quantum number 
i.e. AS'=0 and (c) possible change of the total angular momentum quantum number 
AA=0,±1. 

The selection rules affect the functional shapes of electron impact excitation cross 
sections. However, due to electron exchange effects in electron-impact excitation 
collisions at low energies, these selection rules are not strictly preserved [12]. The most 
signihcant differences of electron impact cross sections are found between singlet-triplet 
(AS'=1) spin-forbidden excitations and singlet-singlet (A5'=0) spin-allowed excitations. 
The excitation cross sections from the ground state (X^S+) to other singlet states 

etc.) are rather insensitive to the impacting electron 
energy in the range >30 eV (Fig. [1]). In the case of spin-forbidden transitions (A5'=l) 
the electron impact excitation occurs via resonance processes. Therefore, the electron 
impact excitation cross sections from the ground state (A^S+) to the triplet states 
(6^S+, c^n„, etc.) are peaked at <20 eV (Fig. H]). 

Altogether the hydrogen molecule has several tens of electronic states [13]. In this 
study singlet states taken into account. This is because the 

electron impact excitation cross sections of and states correspond to over 

75% of the total electron impact excitation cross section to singlet states and cascading 
from the EF^J^g state is the only signihcant effect involving upper states. Other singlet 
states either emit VUV-light in different wavelength range (e.g. and R'^S^) or the 

corresponding electron impact excitation cross sections are orders of magnitude lower 
(e.g. R'lH,, GK^^g) [12]. 

On the other hand, triplet states c^H^, d^Ylu and are taken 

into account in this study. This is due to lack of available cross sections for electron 
impact excitations to other triplet states. The lowest triplet state (&^E+) to which the 
radiative decay chains of and states end, is repulsive. The c^H^ state 
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Figure 1. The cross sections of the most significant electron impact processes of 
hydrogen molecules (y = 0) with hot electrons (i?e > 9 eV). Excitations from the 
ground state to singlet states are marked with black color, to triplet states with red 
color and ionization with blue color. The cross sections are from Ref. [T^ . 


is partially metastable depending on the symmetry (c^n+ or c^H") |T3]. The radiative 
lifetime of the metastable state depends on the vibrational level and ranges from 
1ms (z^=0) to 10 ps (i^=3) [inilTG]. The c^fl^ state is sensitive to electron impact de¬ 
excitation/ionization [T^IIH] and collisional quenching [T9] . 

Plasma emission spectroscopy yields information about population densities of 
molecules on excited states, which are affected by several processes. The diagnostics 
presented in this study are based on the premise, that the VUV-light emission rate is 
equal to the electron impact excitation rate from the ground state to the corresponding 
electronic states (corona model). This sets limitations on the temporal resolution of the 
diagnostic as well as plasma parameters, which are fullhlled by most low temperature 
laboratory plasmas including hydrogen ion sources (see e.g. Refs. mm and references 
therein): 

i) The shortest temporal scale that can be studied is on the order of 1 ns for singlet 
transitions and 100 ns for triplet transitions 

ii) The plasma density is less than 10^^ cm“^ 
hi) The neutral gas pressure is under 500 Pa 

iv) The ion temperature Tj is lower than the electron temperature Tg 

The listed limitations do not concern the metastable c^fl^ state, which is discussed 
thoroughly in Sections 13.31 and 13.41 
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The best possible temporal resolution of the diagnostic method is limited by the 
radiative lifetime of the excited states (delay in spontaneous emission). The radiative 
lifetimes of the lowest singlet states and and lowest triplet states are 

under 1 ns and 10-40 ns respectively [T6] . 

The plasma density is limited by collisional de-excitation and ionization rate from 
the excited states. Signihcant collisional de-excitation of state requires the 

electron density to be lO^^-lO^® cm“^ [21]. For the state the corresponding number 
is 10^^ cm' Since the cross sections of electron-impact ionization from the excited 

states and are on the same order of magnitude (as a maximum) with the 

electron impact de-excitation cross sections iia. the given maximum plasma densities 
are valid also for ionization. 

Electronic states can also be de-excited e.g. via excitation transfer with neutral 
particles, charge exchange with ions and penning ionization. The cross sections of these 
processes are, however, known poorly [12]. Most of these de-excitation processes are 
resonant by nature and hence their cross sections could be large, but the required 
collision energies are on the order of a few eV. Thus, they are negligible in low 
temperature plasmas, where the the average energies of neutral particles and ions are 
low. The excitation transfer with neutral particles could limit the maximum neutral gas 
pressure. However, there are no cross section data available for excitation transfer of 
hydrogen molecule. The lower limit for the maximum neutral gas pressure (r=300 K) 
can be estimated by using the excitation transfer cross section of hydrogen atom 
(approximately 10“^^ cm^), which yields a pressure of 5 mbar assuming 10 ns lifetime 
for the excited state. The given cross section is probably several orders of magnitude 
too large (and calculated upper limit of the pressure too low), because of the vibrational 
distribution of the molecules and the Franck-Condon principle decreasing the excitation 
transfer probability of the molecule in comparison to the atom. 

The applied model requires that the electron impact excitation is the dominant 
excitation process. This in turn reqnires that the ion temperature is lower than the 
electron temperature. The effects of photonic excitation and recombination on the 
population densities of excited molecules are negligible in comparison to corresponding 
effects in the case of atomic emission. Molecular hydrogen (laboratory) plasmas are 
optically thin for molecular emission, because the vibrational distribution of molecules 
and Franck-Condon principle result to low probability of photon excitation. On the 
other hand, recombination of the hydrogen molecule leads mainly to dissociation of the 
molecule [12], i.e. 

e + {H ;*-^ H{ls) + H{n > 2), (1) 

which does not contribute to the population densities of excited molecules. 

3. Estimating rate coefficients from VUV-emission 

A typical VUV spectrum of a hydrogen plasma is presented in Fig. |2j The spectrum 
can be divided into four different regions based on the origin of the radiation. Lyman- 
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Figure 2. VUV emission spectrum of a hydrogen plasma of an arc discharge with 840 
W discharge power and 0.5 Pa pressure. The spectrum is not correlated for spectral 
transmittance. 

series (mainly Lyman-alpha) radiation is emitted by hydrogen atom while three 
dominant molecular transitions emit VUV-light. Those are Lyman- and Werner-bands, 
corresponding to transitions from the lowest excited singlet states and C^Tlu) 

to the ground state, and molecular continuum corresponding to the lowest transition 
between triplet states —)-b^S+). Because the lowest triplet state is repulsive, a 

continuum is observed instead of a band structure. 

The Lyman-band and molecular continuum emissions are the most suitable 
transitions for the purpose of plasma diagnostics. This is because signihcant parts 
of Lyman-band and molecular continuum radiation are emitted in ranges of 145- 
165 nm and 170-240 nm including negligible contribution from other transitions. 
Furthermore, structure of Lyman-band and molecular continuum emissions do not 
depend (signihcantly) on the plasma parameters and cascade effects from the upper 
states to state are well known, which allows straightforward estimation of total 

emission from the measured signals in the specihc ranges. 

VUV-emission can be used as a probe for the most important plasma processes, 
such as ionization rate, vibrational excitation rate to high vibrational levels, production 
rate of metastable c^Hu states, molecular dissociation rate and hot electron density. 
These rates are sensitive to distributions of electron energy and molecular vibrational 
level. 

The vibrational distribution of neutral hydrogen molecules is not in the thermal 
equilibrium, especially at high vibrational levels. However, a majority (>90%) of neutral 
molecules in typical laboratory plasmas are on the lowest vibrational levels {v < 4) (e.g. 
Refs. [22l|23] and references therein). This part of the distribution can be described 
reasonably well with the Boltzmann distribution, i.e. using a vibrational temperature 
Tyib. The vibrational temperature of the lowest vibrational levels has been found to be 
between 100 K and 10000 K (e.g. Refs. [23l|25 and references therein). These numbers 
are used as extremes when estimating the uncertainty of the presented diagnostic caused 
by the vibrational distribution. This is because the cross sections of electron impact 
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ionization and excitation to electronic states are less than linearly proportional to 
the vibrational quantum number. Molecules have also rotational levels, which have 
a negligible effect on the discussed electron impact processes due to minimal energy 
exchange (~0.01 eV) in transitions and long collision times (frozen rotation) [12]. 

The electron energy distribution functions (EEDF) of low temperature plasmas are 
typically non-Maxwellian. The distributions are often bi-Maxwellian at low pressures 
and Druyvesteyn-like at high pressures [25]. At high plasma density (and ionization 
degree) electron-electron collisions tend to drive the EEDF towards a Maxwellian [25] . 
Furthermore, superpositions of Maxwellian and flat (see definition later) distributions 
have been commonly found in filament arc discharges (e.g. Refs. [261 - 128] ). When 
electron-electron collisions dominate in the entire energy range, the resulting EEDF is 
Maxwellian by definition. In hydrogen plasmas this is the situation only at low electron 
energies (Eg<9 eV). The EEDF at high electron energies is determined by elastic and 
inelastic collisions, plasma heating method and electron confinement. 

Typically only high energy electrons (i7e>20 eV) contribute to ionization and 
excitation to singlet states. For example, in the case of the EEDFs reported in 
Refs. |29H3T], >20 eV electrons contribute more than 99% to the ionization and 76- 
91% to the excitation to state although their density is only 1-3% of the total 

electron density in Ref. [2^130] or 30% in Ref. [3T]. Therefore, ionization and singlet 
state excitation rates can be described reasonably accurately by focusing only to the 
high energy part of the EEDF. Excitations to triplet states are most sensitive to the 
electron energy in the range of 10-20 eV. 

In this paper two EEDFs have been chosen for a closer study from a mathematical 
point of view to describe the distribution of electrons at energies exceeding the excitation 
and ionization thresholds (Fig. [3]). At high energies the Maxwellian distribution and the 
tail of the bi-Maxwellian distribution can be described by the Boltzmann probability 
function 


f{Ee) = 


( 2 ) 


where A is a normalization factor, E^. is electron energy, is the Boltzmann constant 
and Tf. is electron temperature. The high-energy tail of the EEDF in filament arc 
discharges can be described by a flat distribution of energies. 



( 3 ) 


where A is a normalization factor and E^ax is the end point (or maximum) energy 


of the plasma electrons corresponding to the potential difference between the cathode 


and the plasma. The rate coefficients of these two distributions differ significantly. If 
the average energy of the distribution {^kbTe for Maxwellian distribution and E^ax/^^ 
for the flat distribution) changes, the electron density at energies corresponding to the 
most sensitive range of the rate coefficients (6 eV<i7e<40 eV, see discussion in Sections 
13.1113.411 changes exponentially in the case of Maxwellian distribution and lineraly in the 
case of the flat distribution. 




8 



Figure 3. Examples of electron energy density functions used in this study. Hot 
electron versions of distributions are used for calculating the rate coefficients presented 
in Section [3^ by renormalizing the Ef, > 9 eV part. The choice of 9 eV corresponds 
to the threshold energy of —>• electron impact excitation for vibrational 

levels 1 / < 4. 


The VUV-emission rate corresponds to the excitation rate of neutral molecules, if 
the plasma parameters are within the limits specihed in section [2l Comparison of the 
rate coefficients of different processes allows estimating their volumetric reaction rates 
from the measured VUV-emission rates of specihc emission bands. The volumetric rate 
R of an electron impact process can be described as 


R = rieUn 


f{v)va{v)dv = HnUe {v(t) , 


( 4 ) 


where Ug is the electron density, the neutral molecule density, f normalized electron 
velocity distribution function (EVDF), v the electron velocity and a the cross section 
of the process. The term (va) is the rate coefficient which depends on the functional 
shapes of the EVDF and the cross section. 

In molecular plasmas the volumetric rate of any given process depends on the 
vibrational distribution. The total volumetric rate can be calculated as a weighted 
average of the rate coefficients corresponding to individual vibrational levels. For 
Boltzmann distribution the result can be written as 


R(d^vib) '^n^i 


^.(w).exp(- 


kT,,, 


Eiexp(- 




( 5 ) 


kTvib ^ 


where is the vibrational temperature, Ei is the energy of the vibrational level i on 
the ground state and is the total excitation rate coefficient from the ground state 
vibrational level i to all possible vibrational levels of the upper state. 

If the volumetric rate of a specihc process (i?i) can be determined e.g. by measuring 
the volumetric emission rate of a specihc band of VUV-light, the volumetric rate of 
another process {R 2 ) can be estimated from 


Ti^?7.eCV2 


( 6 ) 
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If the same neutral species (n„) is involved in both processes, the electron and neutral 
densities cancel out. Thus, the volumetric rate R 2 can be obtained by multiplying the 
measured volumetric rate Ri with the corresponding rate coefficients, i.e. 

(Tuil) 


R 2 = 


-Ri — K2^i{Tyib)Ri, 


( 7 ) 


where X 2 ,i is the ratio of the rate coefficients. 

It is of note that the ratio of the rate coefficients, is signihcantly less sensitive to 
changes of the EEDF (e.g. in terms of average electron energy or functional shape), and 
often to the vibrational temperature, than individual rate coefficients. The sensitivity 
of the ii' 2,1 to the EEDF can be understood by studying the cross sections of the 
corresponding processes (Fig. 1) and their ratios at different electron energies. The 
rate coefficient ratio is sensitive to changes of the EEDF in the energy range where 
both, the ratio of the cross sections changes signihcantly, and the absolute values of the 
cross sections are large. Typically the sensitive energy range is a narrow band close to 
the threshold energies. For example, when the excitation rate coefficients of similar types 
of transitions (in terms of the selection rules) are compared the sensitive energy range 
is 10-20 eV. In such cases (similar types of transitions) the vibrational temperature 
dependence of ^ 2,1 is reduced due to similar dependence of the individual cross sections 
on it. Therefore, in some cases ^ 2,1 is insensitive to the average electron energy, shape 
of EEDF and vibrational temperature although individual rate coefficients can vary 
by orders of magnitude within the same ranges of EEDF and vibrational temperature 
variations. Such situations are presented in Sections 3.2-3.6, where the extreme values 
of K 2 ^i vary only by <30% in comparison to their average values in the entire energy 
range of (both) EEDFs and vibrational temperatures. 

There are several reviews about electron impact cross sections in low temperature 
hydrogen plasmas (e.g. Refs. [T2ll32ll35] ). Data from Ref. [12] has been used in this study 
because it includes the most complete critically assessed set of cross sections including 
htting functions in analytic form. 


3.1. Molecular ionization rate 

Hydrogen molecules can be ionized by electron impact non-dissociatively via 
state and dissociatively via i/ 2 ^(X^S+; e') and H 2 e) states [T2] . 

The cross-section of non-dissociative ionization is more than an order of magnitude 
larger {v = 0) than the total cross section of the dissociative ionization [12]. Hence, 
only non-dissociative ionization is taken into account in this study. The functional 
shapes of electron impact cross sections of ionization and spin-allowed excitations to the 
electronic states are similar. The molecular ionization rate can be estimated from the 
measured Lyman band emission rate with Eq. [ 7 ] by estimating the ratio of the ionization 
and —)■ excitation rate coefficients, Kion,B, plotted for different EEDFs and 

vibrational temperatures in Figs. |4]and[5l 

In the case of Kion,B, the most sensitive range is 10-40 eV. The sensitivity in 
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Figure 4. The ratio Kion,B of ionization and ^ excitation rate 

coefficients as a function of the temperature (Tg) of a Maxwellian EEDF and vibrational 
temperature of the hydrogen molecules. Calculated using cross sections from Ref. [12]. 



Figure 5. The ratio Kion,B of ionization and —)> excitation rate 

coefficients as a function of the maximum energy {Emax) of a flat EEDF and vibrational 
temperature of the hydrogen molecules. Calculated using cross sections from Ref. m- 


this range is mostly due to different threshold energies of the processes; 11.6 eV for 
i?^S+-excitation and 15.4 eV for ionization. In the case of Maxwellian EEDF, there 
is a significant plasma temperature dependence of the rate coefficient ratio due to 
exponential nature of the EEDF in the range of 10-40 eV. In the case of the flat EEDF, 
the ratio of the rate coefficients does not depend strongly on the maximum energy of 
the distribution if Emax>50 eV, which is typical in filament driven arc discharges. 

It can be concluded from Figs. [Hand (5] that the uncertainty of ii"ion,s caused by 
the vibrational temperature is <±15%, if T^ift=6000 K is used and the real vibrational 
temperature is in the range of 100-10000 K. If the electron temperature can be estimated 
(or measured) with an accuracy of ±20% which can be considered as a typical number 
[3T| . it causes an uncertainty of <±20% to Kion,B- Hence, the total uncertainty of 
Kion,B is <±25% in the case of Maxwellian EEDF. It must be emphasized that the flat 
distribution does not accurately describe the real EEDF. However, the uncertainty of 
the diagnostics results caused by such deviation can be described as a variation of (the 
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Figure 6. The ratio Kbc,b of C^IIu and —>■ excitation 

rate coefficients as a function of the temperature (Tg) of a Maxwellian EEDF and 
vibrational temperature of the hydrogen molecules. Calculated using cross sections 
from Ref. |12) . 


effective) Emax of an ideal flat EEDF. This can be demonstrated as follows: (1) choosing 
Kion,B = 2 corresponds to an nncertainty of <±25% for 60 eV<£'maa:<250 eV or (2) 
choosing Kion,B = 1-5 corresponds to an nncertainty of <±35% for 35 eV<Ejnax<85 eV. 

3.2. Excitation rate to and C^Tlu states 

The vibrational distribntion of low temperatnre molecnlar hydrogen plasmas is not 
in thermal eqnilibrinm (e.g. Refs. |22l[23] and references therein). The fraction of 
molecnles at high vibrational levels affects the rates of several plasma processes, e.g. 
the volnme prodnction of negative hydrogen ions throngh dissociative attachment [1]. 
It has been conclnded [Il[36] , that most of the vibrational excitations popnlating states 
u > 5 are preceded by electron impact excitation to the and C^Hu singlet 

states. Radiative transitions from the given excited electronic states to the gronnd state 
populate high vibrational levels, when the vibrational level changes in the electronic 
transition according to the Franck-Condon principle [37]. By assuming that transitions 
between vibrational levels in electron impact excitation to electronic states follow the 
Franck-Condon factors it can be calculated that 63-67% of the (de-)excitations emitting 
in Lyman-band and 46-49% emitting in Werner-band to ground state vibrational levels 
n > 5. Approximately 15% of transitions from and states lead to such high 

vibrational levels (vibrational continuum) that the molecule dissociates instantly [38] . 
The corresponding dissociation rate does not depend strongly on the electron energy or 
vibrational temperature of the molecules. 

The total excitation rate to and C'^fl^ states can be estimated from the 

measured Lyman-band emission with Eq. [7| by comparing the total excitation rate 
coefficients to B^Y3^ and C'^n^ states to the excitation rate coefficient to B^Y3^ state. 
The ratio of the rate coefficients, Kbc,b, is plotted for different vibrational temperatures 
and EEDFs in Figs. |6]and[71 
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Figure 7. The ratio Kbc,b of C^IIu and —>• i?^E+ excitation 

rate coefficients as a function of the maximum energy {Emax) of a flat EEDF and 
vibrational temperature of the hydrogen molecules. Calculated using cross sections 
from Ref. m- 


The rate coefficient ratio is only weakly sensitive to variations of the EEDF in the 
range of 10-20 eV. This can be explained by the similarity of the functional shapes 
of the cross sections and threshold energies (Fig. [1]). For example, using median 
values of Kbc,b = 1-64 for Maxwellian EEDF or Kbc,b = 1-79 for the flat EEDF 
cause corresponding uncertainties of <±10%, if the real electron temperature is in the 
range of 3 eV-25 eV or the effective Emax is in the range of 30 eV-250 eV. On the 
other hand, using = 6000 K causes an uncertainty of <±4% to Kbc,b if the real 

vibrational temperature is in the range of 100-10000 K. Therefore, the total uncertainty 
corresponding to the selected median values of Kbc,b (in respective ranges of EEDF 
vibrational temperature variations) is less than 11%. 

3.3. Maximum production rate of metastable hydrogen molecules 

The metastable state (c^If^) of the hydrogen molecule has an important role in molecular 
hydrogen plasmas. Metastable molecules have lower ionization potential and different 
electronic transition probabilities than ground state molecules. Therefore, the existence 
of metastable molecules allows ionization and de-excitation by lower energy electrons, 
which in turn complicates plasma diagnostics based on triplet transitions. Metastable 
states can also affect the production of negative hydrogen ions [1]. 

From the diagnostics point-of-view there are signihcant uncertainties related to the 
metastable c^Iiu state. Different symmetries of the c^If^ state (c^n+ and c^n~) have 
signihcantly different lifetimes. The lifetime of the state is 6.2 ns i.e. only the 

c^II” state can be considered metastable [16]. There are only few studies about the 
state with a signihcant discrepancy between calculations and experiments [T6] . 
The state is argued to be metastable without any distinction between the parities 
of the wavefunction (c^n~) in the electron impact cross-section data (Ref. [12] and 
references therein). Because there are no cross section data about electron impact 
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Figure 8. Ratios of rate coefficients of total triplet (a^E+, c^n„, 

excitation to excitation leading to molecular continuum emission {Ktripiet,Mc) and 
metastable state excitation to excitation leading to molecular continuum emission 
as a function of temperature (Te) of a Maxwellian EEDF. The given range 
corresponds to the assumption that either all (min) or none (max) excitations to the 
metastable c^n„ state lead to emission within the molecular continuum. Calculated 
using cross sections from Ref. [T^ . 


excitation to the state, the state is assumed to behave as the metastable 

c^II” state. Moreover, experiments and theoretical calculations about electron impact 
excitation cross sections from the ground state to c^If^ state differ by a factor 

of 2-3 [l2]. There is also signihcant difference between calculated cross sections e.g. in 
Refs. [12] and [39] . 

From the point-of-view of molecular continuum emission there are three different 
types of triplet states. The lowest triplet state, is repulsive, and electron impact 

excitation to the this state leads to non-radiative dissociation. Triplet states d^Hu 

and lead to emission within the molecular continuum. The metastable state c^IIu 
can either spontaneously decay, transform to radiative triplet state (a^S+) by electron 
impact de-excitation or collisional quenching, the molecule can drift to the wall of the 
plasma chamber or the molecule can be ionized by electron impact. Therefore, the 
contribution of the metastable molecules on the molecular continuum emission depends 
on the plasma parameters and the dimensions of the plasma chamber. 

The volumetric production rate of metastable molecules (c^II”) can be estimated 
from the measured volumetric emission rate of molecular continuum —)■ fe^S+) 

radiation and the ratio of the rate coefficients of excitation to metastable states and 
total excitation to triplet states leading to molecular continuum emission. The minimum 
production rate of metastable molecules can be estimated if most of the excitations to 
c^n.u-state are assumed to eventually emit a photon within the molecular continuum 
(e.g. via collisional quenching or decay to a^S+). The maximum of metastable state 
production rate corresponds to the situation in which all metastable states decay 
without photon emission in the molecular continuum. Because of the large cross 
section of electron impact excitation from the ground state to the state, there 
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Figure 9. Ratios of rate coefficients of total triplet excitation (a^E+, c^IIu, 

e^E+) to excitation leading to molecular continuum emission (Ktripiet,Mc) and 
metastable state excitation to excitation leading to molecular continuum emission 
as a function of maximum energy (Emax) of a flat EEDF. The given range 
corresponds to the assumption that either all (min) or none (max) excitations to the 
metastable c^n„ state lead to emission within the molecular continuum. Calculated 
using cross sections from Ref. [T^ . 


is approximately a factor of two difference between these two extremes. 

The nncertainty cansed by snch approximations (extreme cases) can be reduced 
by comparing the transition probabilities of different processes affecting the population 
density of the c^If^ state. Typically, the most signihcant de-excitation processes are 
electron impact de-excitation, collisional quenching with neutral particles and collision 
to the wall of the plasma chamber, as discussed in Ref. [ 10 ] • The radiative lifetime of the 
c^II” state depends on the vibrational level and ranges from 1 ms {i^=0) to 10 ps (z^=3) 
[41] . These radiative lifetimes correspond to de-excitation transition probabilities of 
10^ 1/s (z^=0) and 10® 1/s (z^=3), respectively. The collisional quenching rate coefficient 
of the c^II” state in collisions with hydrogen molecules at 300 K is 1.88 x 10“® cm®/s [T9] . 
This means that the transition probability of collisional quenching depends on the 
neutral gas pressure as 4.7 x 10®l/sPa. The rate coefficient of the most significant 
electron impact de-excitation process (the superelastic transition c^Ylu —t b^^u) is on 
the order of 10“^cm^/s regardless of the electron temperature HZ]. Thus, collisional 
de-excitation is significant if the neutral gas pressure is more than 2 x 10“^ Pa and/or 
the electron density is more than 10^°cm“® for c^n„(i^=0) or 0.2 Pa and/or 10^® cm“^ 
for c®n„(n = 3) state. Furthermore, the c^If” state de-excites due to wall collisions [TO] . 
This is a signihcant process in small plasma sources at low neutral gas pressure and 
electron density, because the average thermal velocity of hydrogen molecules (T=300 K, 
1.9mmps“^) corresponds to a distance of some centimetres during the radiative lifetime 
of the c^II” state. 

The ratio of the total rate coefficient of electron impact excitation to metastable 
state (c^IIu) and the total rate coefficient of electron impact excitation to the molecular 
continuum radiative states (a^S+, (PHu, with/without c^n^), and 
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^c^n,MCmax) plotted for different EEDFs in Figs. [H] and IHl The insensitivity of the 
rate coefficient ratio to the EEDF can be explained by similar functional shapes and 
threshold energies of the processes. The contribution of different vibrational levels is 
not taken into account due to the lack of numerical data [T^[39] . However, it can be 
argued that because increasing vibrational level decreases the threshold energies and 
increases the cross sections almost identically for electron impact excitation from the 
ground state to any triplet state [39], the effect of vibrational temperature on 
is almost negligible (of the same order as in Figs. [6] and [7]). 

Using median values of A'c3n,MCmax=0-60/i^c3n,Mc^i„=l-55 (Maxwellian EEDF) 
and Kc3n,MCmax=0-56/7Cc3n,MC™i„=l-26 (flat EEDF) cause an uncertainty of <17% if 
the electron temperature is in the range of 3-25 eV or <5% if the effective Emax is in the 
range of 30-250 eV. The most significant uncertainties (that depend on each other) are, 
however, caused by the contribution of metastable molecules (c^n+) to the molecular 
continuum emission and the uncertainty of the cross section data. In the worst case, 
both of them can cause a factor of two deviation between the diagnostics result and the 
reality. The most signihcant factor is the uncertainty of the ratio of excitation cross 
section to metastable state 77^ —)■ c^H^ and total excitation cross section to other 
triplet states —)■ (a^S+, d^H^, e^S+). 

3.4- Molecule dissociation rate via triplet state excitation 

There are several processes which lead to molecule dissociation in hydrogen plasmas 
e.g. dissociative electron attachment, recombination, ionization, excitation to triplet 
states, excitation to singlet states on high vibrational level and processes of molecular 
hydrogen ion [T2|. The significance of each process depends on the plasma parameters. 
Unfortunately, there are no reliable experimental data for the cross sections of molecule 
dissociation as discussed in Ref [33] . It has been argued that electron impact excitation 
to triplet states is the main dissociative channel for the hydrogen molecule [331IM] . 

The minimum dissociation rate of the molecules by the lowest, repulsive, triplet 
state can be estimated from the measured molecular continuum emission. In this 

case, the rate coefficient of electron impact excitation to all triplet states is compared to 
the rate coefficient of electron impact excitation to the states emitting a photon within 
the molecular continuum. It can be argued that all of the electron impact excitations 
to the state lead to the dissociation of the molecule, because both, the radiative 
and collisional de-excitation of c^H^i lead mainly to state. 

Two extremes of the rate coefficient ratios can be considered by assuming that 
either all {Ktripiet,MC„,iJ none {Ktripiet,MCmaJ of the metastable states lead to 
molecular continuum emission. There is approximately a factor of two difference between 
these extremes. The coefficient ratios Ktripiet,MCmin and Ktriplet,MCma^ are plotted for 
Maxwellian and flat EEDFs in Figs. [8]and[9l The error analysis discussed in Section 
13.31 for the production rate of metastable c^H^i is valid also for dissociation rate of 
hydrogen molecules discussed in this section. 
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Figure 10. The rate coefficient {av) of electron impact excitation 

for a Maxwellian EEDF as a function of electron temperature (Tg). Calculated using 
cross sections from Ref. |12) . 


3.5. Hot electron density 


The density of hot electrons {Eg > 9 eV) can be negligible in comparison to the total 
electron density of the plasma. However, it is the only part of the EEDF which prodnces 
(positive) ions and electrons from the nentral gas. Therefore, the diagnostic of this 
population is important especially from the point-of-view of applications, e.g. optimizing 
the performance of ion sources. The {av) of —)■ electron impact excitation 

is rather insensitive to the electron energy, which allows estimating the hot electron 
density. When the density of neutral molecules is known and the ionization degree 
is low enough not to affect the neutral density, the hot electron density can estimated 
with 


B'T.q, - 


n. = 


Ly—band 


Ur 


(va) 


( 8 ) 


where Riy-band is the measured (volumetric) emission rate of Lyman-band radiation 
and {va) is the rate coefficient. The rate coefficient for Maxwellian and flat EEDFs is 
presented in Figs. fTOlandfTTl 

It can be concluded from Figs. [Tn]and[TTl that the uncertainty of {va) caused by 
the vibrational temperature varies from 40% (Te=3 eV) to 20% (Te=24 eV) in the case 
of Maxwellian distribution and from 21% [Emax=^^ eV) to 13% (i?maa;=250 eV) in the 
case of flat distribution if Tyn, = 6000 K is used and the real vibrational temperature is in 
the range of 100-10000 K. If the electron temperature can be estimated (or measured) 
with an accuracy of ±20%, it causes an uncertainty of <±15% to {va). Using the 
median value of {va) =1.6 x 10“® cm^/s causes an uncertainty of <±33% if the effective 
Ejnax is in the range of 30-250 eV. Thus, the total uncertainty of {va) is <±43% for 
Maxwellian EEDF and <±39% for the flat EEDF in the described situations. 
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Figure 11. The rate coefficient {av) of ^ electron impact excitation for 

a flat EEDF as a function of the maximum electron energy (Emax)- Calculated using 
cross sections from Ref. m- 


3.6. Error analysis 

The presented methods to analyse reaction rates in molecular hydrogen plasma are based 
on absolute measurements of the volumetric emission rates of specific emission bands. 
There are several sources of error affecting the estimated absolute reaction rates: 

i) measurement technique 

ii) overlapping emission bands 

iii) coverage of emission band by the measurement 

iv) cascade from upper states 

v) deviation of the assumed EEDF from the real one 

vi) estimation of the vibrational distribution 

vii) effect of other plasma processes e.g. collisions of metastable molecules 

viii) uncertainty of the cross sections 

ix) validity of the assumptions listed in section |2] 

In an earlier paper m, it was estimated that the uncertainty of the measured 
volumetric emission rates of specific emission band in specihc wavelength range is less 
than 18% (the experimental setup is described briefly in section 0]). This error covers 
both, the error of the measurement itself (item H]) and the overlap of the emission bands 
(item m]) in the selected wavelength ranges of 145-170 nm for Lyman-band and 170- 
240 nm for molecular continuum. 

The measured ranges of Lyman-band and molecular continuum emission do not 
cover the entire emission band. The coverage can be estimated by studying synthetic 
emission spectra. The synthetic spectrum can be calculated, if vibrational the 
distribution of the upper state and the Franck-Condon factors governing the change 
of vibrational level in electronic transitions are known. The vibrational distribution 
of the state can be estimated by using distributions described in Ref. w or 
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by assuming that the vibrational temperature of ground state hydrogen molecules is 
between 0-10000 K and that vibrational excitation in the electron impact excitations to 
the electronic state follows the Franck-Condon factors (hrst order approximation). 
Using the afore-mentioned vibrational distributions at state and Franck-Condon 

factors from Ref. 133 yields that the range of 145-170 nm covers 38-43% of the total 
Lyman-band emission. Using the synthetic spectrum of molecular continuum from 
Ref. [23], yields that the range 170-240 nm covers 59-62% of the total molecular 
continuum emission if the vibrational temperature of the neutral gas is in the range 
of 0-9000 K. If median values (40% and 61% respectively) are used, the variation of 
the correction factors correspond to uncertainties of 5% and 3% for Lyman-band and 
molecular continuum emissions, respectively. 

Recent cross section data of electron impact excitation to EF^'Eg state |43| imply 
that 9-15% of the total Lyman-band emission is caused by cascade from the EF^Tig 
state in the case of the studied EEDFs. Synthetic spectra =100-10000 K) imply, 
that the measured signal (145-170 nm) includes 31-34% of the total cascade effect. 
Thus, 86-93% of the measured signal is caused by direct excitation to state. 

Using a median value (89.5%) for the correction includes a relative uncertainty of 4% 
(estimated similar to the previous paragraph). Cascading from upper states affects also 
the molecular continuum. Without accurate cross sections the correction due to cascade 
contribution (from states which are not taken into account in this study) to molecular 
continuum emission can not be estimated accurately. However, it can be argued that 
the cascade effect is small (i.e. <10%), because the cross sections of electron impact 
decrease signihcantly as a function of the main quantum number. 

The total uncertainty caused by items IMvl could be reduced at least by a factor 
of two with improved calibration of the photodiode and more detailed analysis of 
overlapping emission bands, hlter transmission and stray emission. 

The uncertainty caused by estimating the EEDF and vibrational temperature 
(items |y] and |yi|) are discussed in Sections 13.1113.51 The uncertainty is ll%-25% in 
the case of listed plasma processes and 26-43% in the case of the hot electron density. 
These values correspond to situations in which the vibrational distribution of the plasma 
is unknown and only rough approximation of the EEDF is available. Therefore, detailed 
knowledge on the EEDF and vibrational distribution would reduce these uncertainties 
signihcantly. 

The lack of experimental data on relevant cross sections is problematic [T21I33] . The 
calculations of the cross sections have become more accurate, but experimental studies 
include uncertainties of several tens of percents especially at low electron energies [33] . 
Moreover, there is a variation from several tens of percent up to a factor of 2-3 
between calculations especially in the case of excitation to triplet states [T2l[39]. The 
presented diagnostics of plasma processes is mostly sensitive to relative error between 
the cross sections due to the comparison of rate coefficients. The absolute error of the 
cross sections does not have a signihcant effect on the diagnostics results if the main 
contribution is systematic error, which is similar for each cross section. The uncertainty 
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Figure 12. LIISA H ion source and the irradiance meter used in Ref. m 


caused by the cross section data can not be estimated accurately based on the available 
data, but it can be argued that it is probably comparable to or larger than the total 
uncertainty caused by items 

Deviations from the assumptions made in Section [2] can be considered a minor 
source of error. This is because typical ion source plasmas are within the range of these 
assumptions by a great margin (orders of magnitude) [n[20]. 

4. An example of plasma analysis 

Absolute and relative VUV emission have been measured earlier from the hlament 
driven arc discharge negative ion source, LIISA [TT] (Fig. [12]). The plasma chamber 
(diameter 9 cm, length 31 cm) is made of copper, but it is coated with tantalum due 
to the evaporation of the hlament. The arc discharge power was varied in the range 
of 140-840 W by adjusting both, the arc current between 4-12 A and the arc voltage 
between 35-70 V. The presented values were measured at the plasma chamber pressure 
of 0.35 Pa, which corresponds 5 seem feed rate of hydrogen ga^. The relative VUV- 
emission was not observed to depend signihcantly on the pressure in the range between 
1.2 X 10“^ mbar and 2.0 x 10“^ mbar m- The total electron density of a similar ion 
source [15] has been measured to be always below 10^^ 1/cm^. Thus, the requirements 
for the diagnostics described in the Section [2] are valid in the studied ion source. 

The measurements were performed with a VUV-irradiance meter, which consists 
of a factory calibrated photodiode (IRD-inc SXUV20BNC) and optical bandpass Liters 

f The optimum plasma chamber pressure is presented incorrectly in Ref. The correct plasma 

chamber pressure is 0.35 Pa (measured with Pirani gauge including up to factor 2 uncertainty [44]) 
instead of 0.17 Pa . The gas feed rate of 15 seem was estimated in Ref. m- The value 5 seem is 
experimentally verified. 
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Table 1. Absolute VUV-emission of the filament arc discharge m- Relative error for 
emission power is ±18%. 


Emission band 

(A) 

Emission power 

Photon emission 

Lyman-band /Parc 

160 nm 

40 mW/kW cm^ 

3.3 X lO^M/kWscm^ 

Mol. COnt/Jarc 

195 nm 

0.55mW/Acm^ 

5.4 X 10^^ 1/Ascm^ 


of specific wavelength ranges - Lyman-band (161 nm, FWHM 20 nm) and molecular 
continuum (180 nm, FWHM 40 nm). The device was looking into the plasma along an 
axial line-of-sight through the extraction aperture (diameter 9 mm) of the ion source 
from a distance of approximately 1.5 m. The photodiode (an effective area collimated 
to 3.14 mm^) observed a plasma volume of approximately 20 cm^ out of the 2340 cm^ 
total volume of the plasma chamber. The transmittances of the hlters were measured by 
using a spectrometer, which consists of a monochromator (McPherson Model 234/302) 
and a photomultiplier (FT Enterprices 9406B) and utilizing the ion source as a light 
source. The background signals caused by the stray transmittance and overlapping of 
adjacent emission bands were taken into account, and the corresponding uncertainty due 
to such corrections was included to the total uncertainty. Thus, the presented values 
(Table [I]) correspond to Lyman-band emission and molecular continuum emission in the 
wavelength ranges of 145-170 nm and 170-240 nm, respectively. The average wavelengths 
in Table [Ucorrespond to energy weighted average wavelengths of the (corrected) emission 
in the given wavelength range. 

The conhnement of primary electrons (emitted by the hlament) in the multicusp 
magnetic held of hlament arc discharges is known to be ehective (e.g. Refs. [251146] ). 
The same conclusion has been derived from the VUV-emission data obtained with the 
LIISA ion source. It has been observed that the arc power is efficiently absorbed in 
inelastic collisions between primary electrons and neutrals m- Furthermore, it has been 
observed that Lyman band radiation is linearly proportional to the arc discharge power, 
while the molecular continuum radiation is mainly proportional to the arc discharge 
current in the entire parameter range. Therefore, it is convenient to transform the 
measured irradiance values to average volumetric emission rates in the line-of-sight 
volume and normalize the results to the arc power (Lyman-band) or the arc current 
(molecular continuum) as presented in Table [1] 

The EEDF of the hot electron component in hlament arc discharges is best described 
by the hat distribution discussed in Section [3l Such distribution has been observed 
experimentally and deduced numerically (e.g. Refs. [261 - 128] ). Furthermore, the electron 
energies can not exceed the potential diherence of the hlament and the plasma, i.e. the 
Maxwellian distribution would be a non-physical representation of the high energy tail. 

The estimated volumetric rates and hot electron densities in the line-of-sight volume 
are presented in Table [2l The total emission rate of Lyman-band (direct excitation to 
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Table 2. Calculated reaction rates of different processes and hot electron densities. 
Variation of the rates corresponds to 35 - 70 eV maximum energies of the flat EEDF. 
The rates have been normalized with the arc power or arc current depending on their 
response to these parameters. 


Process 

Vibrational temperature 

100 K 10000 K 


Ionization 

8-15 

7-12 

xl0^®l/kWscm3 

Excitation C^H-u) 

13-14 

12-13 

xl0i®l/kWscm3 

Excitation {v > 5) (i?^E+, C^n„) 

00 

1 

1 

o 

xl0i®l/kWscm3 

Dissociation (B^E+, C^n„) 

1.9- 2.0 

bo 

1 

xl0i®l/kWscm3 

Dissociation (via 6^E+), min 

1.8 


xlO^® 1/Ascm^ 

Dissociation (via &^E+), max 

4.2 


xlO^® 1/Ascm^ 

Eliu states, min 

0.5 


xlO^^ 1/Ascm^ 

Eliu states, max 

1.1 


xlO^® 1/Ascm^ 

Hot electron density 

5.7-8.5 

4.1 - 5.7 

xl0^° 1/kW cm^ 


state) and molecular continuum are calculated by using the equation 

p _ tineas , 

^tot — 7 ^cas) 


(9) 


where Rmeas is the measured photon emission rate (Table [I]), kcov is a correction factor 
taking into account the coverage of the measurement compared to the total emission and 
kcas is a correction factor taking into account the cascade effect. Values of 0.4 (Lyman 
band) and 0.6 (molecular continuum) have been used for kcov and 0.895 (Lyman band) 
and 1 (molecular continuum) for kcas as discussed in Section 13.61 The volumetric rates 
were calculated from equation[7]by using rate coefficient ratios from Figs. [5l [7]and|9l The 
dissociation rate and vibrational excitation rate to V^E+(n > 5) via electronic excitation 
to and C^Hu states are calculated by using the fractions between the measured 

and total rates discussed in Section 13.21 The hot electron densities are estimated with 
Eq. IH]by using the rate coefficient from Fig. [11] and neutral gas density corresponding to 
0.35 Pa pressure at room temperature. Variation of the values in each cell corresponds 
to the variation of applied arc voltages (Emax)- Vibrational temperatures 100 K and 
10000 K have been used to demonstrate the effect of the vibrational temperature. 

The maximum density of metastable molecules can be calculated from the estimated 
production rate and average lifetime of the metastable states. The plasma parameters 
(density and neutral gas pressure) are in the range where, in principle, all of the de¬ 
excitation processes listed in Section 13.31 affect the effective lifetime of the c^If^ state. 
However, the maximum effective lifetime of the state can be estimated to be 

26 ps, which corresponds to the time-of-flight of thermal hydrogen molecule across 
the radius of the plasma chamber (5 cm). Hence, a typical maximum production 
rate of 1.1 x 10^®l/scm^ with 10 A arc current corresponds to a maximum density 
of 4 X 10^^1/cm^ molecules excited to the c^H^ state. This in turn corresponds to 
approximately 0.4% of the neutral gas density at 0.35 Pa at room temperature. 
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Table 3. The average number of a single molecule to undergo a certain process. The 
numbers have been normalized with the arc power or arc current. The calculated 
numbers assume a homogeneous and isotropic plasma emission distribution occupying 
the whole chamber. This assumption overestimates the values by not more than 50% 
as explained in Ref. m 


Process 

Vibrational temperature 

100 K 10000 K 


Ionization 

73 - 141 

61 - 117 

1/kW 

Excitation C^Iiu) 

121 - 129 

115-123 

1/kW 

Excitation {v > 5) (i?^I]+, C^IIu) 

70-74 

67-70 

1/kW 

Dissociation (R^S+, C^IIu) 

18-19 

17-18 

1/kW 

Dissociation (via min 

1.7 


1/A 

Dissociation (via max 

3.9 


1/A 

c^n„ states, min 

0.5 


1/A 

c^n„ states, max 

1.0 


1/A 


The total reaction rates can be compared to the neutral gas density, the gas feed 
rate and the number of electrons emitted by the hlament. This allows estimating the 
average number of different reactions experienced by a single molecule during the time 
it spends in the plasma chamber (Table [3]) and the average number of different types of 
inelastic collisions caused by a single electron emitted by the hlament (Table H]). The 
total reaction rates can be estimated by assuming homogeneous and isotropic plasma 
emission prohle. This overestimates the obtained values by less than 50% [TT] . 

The gas feed rate of 5 seem corresponds to 2.1 x 10^® molecules/s. Comparing 
this number to the total number of molecules in the chamber (0.35 Pa of ideal gas, 
N=2.6 X 10^^ molecules) yields an average passage time of 83 ms through the plasma 
chamber. This means that each molecule undergoes on average several tens of ionization 
and singlet state excitation reactions and dissociates several times (Table [3]). 

It is possible to estimate the energy efficiency of the plasma source by comparing 
the total reaction rates to the heating power and the number of electrons emitted 
by the hlament. This is presented in Table IH Each kilowatt of discharge power 
produces approximately 1.3-3.0 x 10^*^ of electrons and ions per second via ionization 
corresponding to 20-40 Amperes of electrical current. Approximately the same number 
of neutral molecules are excited to singlet states per kilowatt of discharge power. The 
molecular continuum emission is proportional to the arc current, and therefore, it 
is reasonable to compare the molecule dissociation rate and the production rate of 
metastable states to the number of electrons emitted by the hlament. It turns out 
that every electron emitted by the hlament dissociates on average 0.6-1.3 molecules via 
triplet states and excites 0.15-0.4 molecules to the metastable state. The dissociation 
rates via triplet and singlet states are comparable under typical ion source parameters 
(arc voltage 70 V, arc current 12 A), i.e. molecules dissociate only 1.5-3 times more 
frequently via triplet states than singlet states. The results are consistent with the 
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Table 4. The average number of ionization, excitation, dissociation and production 
of metastable states per kilowatt of arc discharge power or electrons emitted by 
the filament. The calculated numbers assume a homogeneous and isotropic plasma 
emission distribution occupying the whole chamber. This assumption overestimates 
the values not more than 50% as explained in Ref. m- 


Process 

Vibrational temperature 

100 K 10000 K 


Ionization 

15-30 

13-24 

xl0i®l/kWs 

Excitation C^Iiu) 

25-27 

24-26 

xl0i®l/kWs 

Excitation {v > 5) C^Ilu) 

15 

14- 15 

xl0i®l/kWs 

Dissociation C^n„) 

3.8-4.0 

3.6-3.8 

xl0i®l/kWs 

Dissociation (via min 

0.6 


events/car-c 

Dissociation (via max 

1.3 


events/car-c 

c^n„ states, min 

0.15 


events/car-c 

c^n„ states, max 

0.35 


events/car-c 


explanation given in Ref. There, the conclnsion was that dnring the thermalization 
process each electron emitted by the hlament passes throngh the optimnm energy 
range for triplet excitation, which explains the molecnlar continnnm radiation being 
proportional to the arc cnrrent. 

5. Discussion 

The principles of robnst and straightforward diagnostic of electron impact processes in 
low temperatnre molecnlar hydrogen plasmas have been presented in this stndy. The 
method could contribute signihcantly to the development of plasma and ion sources 
since the effects of mechanical modihcations, for example, could be connected directly 
to the changes of reaction rates in the plasma. Similar information could be obtained 
by measuring the plasma parameters (neutral gas density, electron density and EEDF) 
and applying collision radiative models. 

The presented methods include signihcantly smaller uncertainty than calculating 
the reaction rates from the measured plasma parameters. This is because applying 
the presented method does not require any information about the electron and neutral 
densities and because the rate coefficient ratio K 2 ^i (Eq. [7]) is less sensitive to the electron 
temperature than the rate coefficient (av) itself. The most signihcant difference in the 
uncertainties of different methods is caused by the electron temperature dependence 
of K 2 ^i and {<Jv). This is demonstrated in Fig. [13] showing the normalized ratio of 
ionization and excitation rate coefficients Kion,B and normalized ionization rate 

coefficient as a function of the electron temperature. Although Kion,B is the 

most sensitive A' 2,1 coefficient (discussed in this study) to the variation of the electron 
temperature, it changes only by an order of magnitude in the range of 3 eV<Tf.<25 eV, 
while varies almost three orders of magnitude. For the other processes studied 
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Figure 13. Normalized ionization rate coefficient and ratio of rate coefficients 

{Kion,BY.) as a function of electron temperature in different vibrational temperatures. 


in this paper the difference is even more significant. The coefficient Kion,B is more 
sensitive to the variation of the vibrational temperature than However, for 

other processes, i^ 2 ,i is less sensitive to vibrational temperature than (av) due to similar 
dependence of the rate coefficients on the vibrational temperature. 

Based on the error analysis presented in Section 3.6. it is justified to claim that the 
main sources of uncertainty for determining the ionization rate and total excitation rate 
to and C^H^ states are uncertainties of measurement (relative uncertainty ±18%), 

the coverage of the measured emission (±5%), cascade effects (±4%) and determining 
the coefficients Kion,B and Kbc,b (<±25% within the limits of EEDF and vibrational 
temperature variations discussed in Sections 3.1 and 2.1). Applying the general law 
of error propagation for the independent error sources, it can be argued that with the 
presented methods it is possible to measure the ionization rate and total excitation 
rate and C^Ilu with an accuracy of 32% or better. This number does not take 

into account the uncertainty of the cross section data. The both effect of metastable 
states and discrepancies in cross section data (e.g. in Ref.s [12] and [39]) do not 
allow making a general statement about the uncertainty of diagnostics results based 
on molecular continuum emission. However, the EEDF, the vibrational distribution 
and the density of metastable molecules are often characteristic properties of the given 
plasma source. Therefore, the listed uncertainties mostly affect the systematic error of 
the VUV diagnostic. It is possible to measure relative changes with significantly higher 
accuracy in comparison to absolute values, which is often more important for practical 
development of plasma sources. 

Several simulation codes have been developed for molecular hydrogen plasmas 
during the last decades, e.g. Ref. |23] and references therein. The presented diagnostics 
of the reaction rates could be a powerful tool to benchmark the simulation codes, as 
introduced in Ref. m- The VUV-emission is a consequence of the most significant 
plasma processes (with the largest cross sections), which have been taken account 
by most of the codes. Therefore it should be possible to extract the volumetric 
VUV emission rates as simulation outputs. Comparison of the simulations and 
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the VUV-spectroscopy could validate input parameters of the simulations, e.g. the 
functional shape of the hot part of the EEDF and/or the vibrational distribution. 
The presented diagnostics could allow developing new type of plasma simulations in 
which the measured reaction rates of hot electron processes (ionization, dissociation 
and vibrational excitation rates) could be used as a parameter and the simulation could 
focus on cold electron processes such as plasma diffusion, recombination and dissociative 
electron attachment. 

The presented methods have been applied to analyze a hlament arc discharge. 
Although the plasma parameters depend on mechanical design of the plasma source, 
some results can be generalized. It is well known, that the conhnement of primary 
electrons emitted by the hlament is good in multicusp arc discharges and their energy 
dissipates mainly in inelastic collisions with neutrals. Therefore, the reaction rates 
normalized with respect to the arc power and the arc current (Table 01) should be 
general properties of this type of plasma sources. 

Excitation rate to vibrational levels {u > 5) via electron impact excitation to 
and state has been estimated earlier by Graham [9] to be on the order 

of 10^^ cm“^s“^kW“^. This is approximately an order of magnitude lower than the 
estimated value given in this paper. The deviation is probably caused by the difference 
in the conhnement of the primary electrons since the apparatus used in Ref. [9] does 
not contain a multicusp magnetic held. The density of metastable c^If” states in a 
multicusp ion source has been measured earlier to be on the order of 10^° cm“^kW“^ |40j . 
This is approximately two orders of magnitude lower than the upper limit estimated in 
this study. The diherence is probably caused by uncertainties related to the estimate 
presented in this study i.e. a) the applied cross section is the total cross section to 
states, b) the contribution of c^II^ state to the molecular continuum emission 
c) the estimated maximum ehective lifetime of the c^II” state and d) the assumption 
of homogeneous and isotropic plasma emission prohle. Thus, the estimated maximum 
c^fl” density presented in this study is consistent with the measured value given in 
Ref. 00] . 

The benehts of VUV-spectroscopy arise from the quantum mechanical properties 
of the hydrogen molecule. Processes, like excitation to triplet states leading to molecule 
dissociation and excitation to singlet states leading to production of high vibrational 
levels on the ground state, are directly related to the properties of the hydrogen molecule 
orbitals and Franck-Condon principle. Although those properties are not valid for 
other elements, there are often interesting (chemical) phenomena related to electronic 
excitations and forbidden transitions, e.g. the production of the metastable helium [48] . 
Therefore, it should be possible to utilize the same principle, namely the comparison of 
the rate coefficients with similar functional shapes, for other plasmas as well. 
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